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Abstract 



O ■ The model for constituent glue is presented starting from the perturbation theory in 

nonperturbative QCD vacuum. Green function is constructed for a system containing qq 
qj \ pair and gluon propagating in the vacuum background field, and Hamiltonian approach 

to the problem is formulated. The masses of lowest qqg hybrids with light quarks are 
evaluated, and implication of results for the spectroscopy of exotic states is discussed. 

There is no doubts now that the gluonic degrees of freedom are not only present in the QCD 
Lagrangian, but should also exhibit themselves at the constituent level in such a way that the 
pure gluonium and hybrid states should exist. Experimental candidates for such exotic states 
appear from time to time, causing quite natural excitement in the situation when the standard 
qq nonets are already overpopulated. On the other hand, the state-of-art of the QCD in the 
strong coupling nonperturbative regime still does not allow to predict unambigiously the mass 
scale of the non-gg states. 

Estimates for the masses of the hybrid qqg mesons were carried out within the various QCD- 
inspired approaches, such as bag models [1], constituent gluon models [2], flux-tube model [3,4], 
and QCD sum rules [5], with rather controversial results. For example, the predictions for the 
masses of lightest hybrids with u and d quarks vary from 1.5 GeV to 2.5 GeV. Having in mind 
that model approaches usually contain both theoretical and parameter uncertainties, we claim 
that more direct QCD-motivated treatment is desirable to handle the problem. 

In the present letter we study the qqg system in the framework of Vacuum Background 
Correlator method [6], constructing the hadronic Green functions starting from the QCD La- 
grangian, with the dynamics of interaction defined by the averages of corresponding Wilson 
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loop operators. It may be shown [6,7] that under reasonable assumptions on the behaviour 
of background vacuum correlators the area law asymptotics for Wilson loop averages at large 
distances appears, giving rise to the string-type interaction of constituents. The qq system 
with this interaction was described in [8], and the picture of hybrids in this approach was first 
considered in [9]. 

We introduce a model for constituent glue starting from the perturbation theory expansion in 
nonperturbative QCD vacuum [10]. As it is usually done in the background field formalism, we 
split the gluonic field into the background part and the perturbation a M over background: 

A^ = B^ + a„. (1) 

Ascribing the inhomogeneous part of gauge transformation to the field we can form gauge 
invariant states with the field a M involved. One-gluon hybrid is represented as 

${x q , Xg, Xg) = ^ a (Xg)^(Xg, X g )a^(x g )^J(x g , X q )l[) S (Xq), (2) 

where a. ..5 are the colour indices in the fundamental representation, spin indices are omitted, 
V>(V0 stands for the quark (antiquark) field, a 1 ^ = (X a )^a a , and parallel transporter $ contains 
only background field: 

^{x,y) = Pexp [ X B« p dz, (3) 

The evolution of the state (2) is described by the Green function 

G(x q x s x g , y q y g y g ) =< ^{y q , y q , y g )^ + (x q , x q , x g ) > B , (4) 

where the averaging over background field configurations is implied. 

The Green function for the field a M propagating in the given background field B^ may be 
written in the Eucledian space-time as [10]: 

G~l = D 2 5, u - gF Vil - £>„£>„ + -D^ = (5) 
= + -D^, 

with the background gauge fixing term G a = (D^a^) a , where Df 1 = d\5 ca + gf cha B\. If the 
background field satisfies the equation of motion D^F^ U = 0, then one has M^ U D V = 0, and the 
Green function (5) may be rewritten as 

G, v = {8, v + D,(^)D X • (£ - 1))(D 2 - 2gF)ll. (6) 

The choice £ = corresponds to the Landau gauge, in which the Green function (6) contains 
explicitly the projector P^x onto transverse states:^ 

P^ = S^-D^)D X (7) 



3 The condition D^F^ — is claimed to be not necessary [11] and may be removed using the Coulomb 
background gauge, which has an additional advantage: in this gauge the ghosts can be decoupled explicitly. 
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The next step is to use the Feynman-Schwinger representation [12] of Green function (4) 
to define the effective action for the qqg system. Here we present the simplified version of the 
model, neglecting spin degrees of freedom and reducing the problem to the effective scalar one; 
this approximation corresponds to the omitting the projector in the gluon Green function (6) 
and colour magnetic interaction (term proportional to gF). In a similar way we omit the spin 
dependence in the quark Green function, assuming 

G q = (D*- m y (8) 

Within this approximation the Feynman-Schwinger representation for the Green function 
(4) takes the form 

G(x q x q x g ; y q y q y g ) = (9) 

roc poo poo p 

J ds J ds J dS J DzDzDZexp(-IC q - /Q -K g )<W >b, 

where 

K, q = m 2 s + i J z 2 (T)dr, /Q = m 2 q s + z 2 (r)dT, JC g = ^ J Z 2 (r)dT, 
boundary conditions read 

z(0) = x g , z(0) = x g , Z(0) = x g , 

z(s) = y q , z(s) = y 9 , Z(S) = y g , 
and all the dependence on the gluonic fields B is contained in the Wilson loop operator 

W = \ a ®r q (y g ,x g )K®rA x <i,yq)®ab,r g (yg,x g ), (10) 

the contours T q , T q and T g run over the trajectories of the quarks and the gluon correspondingly 
(see Fig.), $ afe is the transporter in the adjoint representation. 

The Wilson loop configuration (10) can be disentangled using the relation between ordered 
exponents along the gluon path T g in the adjoint and fundamental representation: 

2®ab,T g (x,y) = tr(\ a <f> rg (x,y)\ b <S>r g (y,x)), (11) 
(path T 9 is directed oppositely to T g ), with the result 

<W > B = \ < (W 1 W 2 - jjrW) > B , (12) 

where W\, W2 and W are the Wilson loops in the fundamental representation along the closed 
contours Ci, C2 and C shown at the Figure. The averaging over background can be done using 
the cluster expansion method [6]. The average of two Wilson loops was calculated in [7], and 
it was demonstrated that, assuming the existence of finite gluonic correlation length T g , the 
generalized area law asymptotics may be obtained. For the Wilson loop average (12) one has 
(when both contours C\ and C 2 are large): 

< W >= ^^exp(-a(S 1 + S 2 )), (13) 
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where a is the string tension in the fundamental representation, and S\ and S2 are the minimal 
areas inside the contours C\ and C 2 . The area law (13) holds for practically all the reason- 
able configurations in the qqg system, apart from the special case of the contours C\ and C 2 
embedded into the same plane, when the expression (13) is replaced by 

< W >= ^^exp{-a{S l - S 2 ) - a ad 'S 2 ), Si > S 2 (14) 

o ad i is the string tension in the adjoint representation. Both regimes (13) and (14) match 
smoothly each other at the average distances between the contours C 1 and C 2 of order of 
correlation length T g . Note that to derive the expression (13) we don't use the standard 
procedure [13] of appealing to the limit N c — > 00 and replacing the adjoint line in (10) by the 
double fundamental one. In actual calculations we, however, assume the regime (13) to be 
valid everywhere in the qqg configuration space, arriving in such a way to the picture of the 
qqg hybrid as a system of a gluon with two minimal strings attached, each having a quark (or 
antiquark) at the end. 

To reduce the four-dimensional dynamics in (9) to the three-dimensional one we follow the 
procedure outlined in [8], using the parametrization for which 

z» = (r, f q ), % = (r, fg), = (r, f g ), 

and introducing new dynamical variables 

T T ._ T ■ 

MV = ^o(r), /i 2 (r) = 2=z (t), ^(t) = 

with < r < T. The representation (9) for the Green function is now rewritten as 

G = J DfgDfgDfgDfuD^D^expi-A), (15) 



with the effective action 



+0-^ dpi^Jwfw? - (wi^i) 2 + a J d(3 2 \l wlw' 2 2 - (w 2 w' 2 ) 2 }, 
where the surfaces Si and S 2 are parametrized by the coordinated w^, and 

Assuming the straight-line approximation for the minimal surfaces, 

«V(r,/3i) = z^(t)Pi + Z»(l - Pi), w 2 ^r,f3 2 ) = z^t)(3 2 + Z»(l - ft), 
we arrive to the Lagrangian 

m\ m\ inr\ ji 2 r\ ^f 2 g 
L ~2fii2^ 2 + 2 + 2 + 2 + [<) 
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+<rpi £ dp iy Jl + ll + ap 2 £ dfoyjl + ll 



l\ = {frf q + (1 - f3i)'r g f - 4(A(rVi) + (1 - A)(^Pi)) 2 , 

Pi 

\\ = (p 2 ? q + (i - p 2 y g ) 2 - - 2 m^p2) + (i - fidfafr)) 2 , 

P2 

Pi = r q -r g , p 2 =rq-r g . 

The Lagrangian (17) is the straightforward generalization of the Lagrangian obtained in [8] for 
the qq string. 

To deal with the square root terms in (17) the auxiliary field approach was used in [8]. Here 
we are interested only in the masses of lowest states, so we adopt another strategy, which allows 
simple variational treatment of the problem. 

First, the experience borrowed from the results of papers [8] tells us that for low values of 
relative angular momenta the square roots in (17) can be expanded up to the second order 
in the transverse velocities k, and it was proved to be a good approximation even for the 
massless constituents. This approximation corresponds to the potential-like regime with the 
linear potential V = ap\ + crp 2 , while the terms ~ If (string corrections) can be taken into 
account as a perturbation. With this simplification the Lagrangian (17) becomes quadratic in 
velocities, the centre-of-mass motion is decoupled, and the Hamiltonian (in the Minkowsky 
space-time) can be easily obtained from (17): 

m\ m\ p 1 + p 2 + p 3 p 2 Q 2 
2p 1 2p 2 2 2 n p 2p Q 

where the Jacobi coordinates r and p and conjugated momenta p and Q are introduced: 

r = r q -r q , p = r g ■ 

Pl+ P2 
PlP2 Pz{Pl+P2) 

Pp = — ; — > pq - — ; ; — 

Pi + P2 Pl+ P2 + P3 

- , P2 -> Pi 

Pi = ~P + : r, p2 = ~P ; r. 

Pi +P2 Pi+ P2 

The Hamiltonian (18) contains the fields Pi(r), and one is to integrate over these fields in the 
path integral representation. Instead we proceed in a way suggested in [14]: we find eigenvalues 
£o(Pi) of the Hamiltonian H variationally, assuming pi independent of r, and then minimize 
Eo(Pi) with respect to pi. Such procedure works nicely for the ground states, with the accuracy 
about several percent, and, apart from being very simple technically, allows for the approximate 
solution of the problem of defining the physical transverse states. 

Indeed, we have neglected the spin dependence in the gluon Green function (6). Thus we 
do not expect the eigenf unctions of the Hamiltonian (18) to be transverse and they are really 
not. On the other hand, within the framework of the variational procedure described above we 
can impose the constraint 

^3^0 - P3(f g $) = 0, (19) 
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projecting out the physical hybrid state ty\ = (^q, ( we restore for the moment the gluon spin 
index A in the state (2)). The constraint (19) is compatible with the projector P (U ^(eq.(7))(after 
averaging over background and introducing the variables /ij). In accordance with the condition 
(19) and neglecting string corrections we choose the physical states to be transverse with respect 
to the (three-dimentional) gluon momentum: 

p 3 $ = 0, ^o = 0, (20) 
so that the states contain electric or magnetic gluons : 

^ ~ Y jjm (Q) (21) 



n ~ \j^fn-UQ) + \j ^T[Y jj+lm (Q) (22) 

where j is the total momentum in the gluonic subsystem. (See the paper [15], where the 
gluonium spectrum was calculated within the similar assumptions). 

We stress that the choice (20), though looking quite natural, is nothing but the variational 
ansatz, and, in principle, cannot substitute for rigorous treatment of the problem of transverse 
and longitudinal gluonic degrees of freedom. The latter should be done with inclusion of spin 
variables into the path integral representation for the gluon Green function. 

As the illustrative example of our approach we present the results for the masses of lowest 
hybrids with u and d quarks, assuming the latters to be massless. The radial part of the 
three-body wave function was chosen to be of the simple cluster form ~ exp(— apf — ap2)( u P 
to the centrifugial barrier), and the angular momenta were taken to be the lowest compatible 
with the ansatz (21) or (22). So the ground state hybrid contains a 1 + or Y~ gluon and a 
quark-ant iquark pair with relative angular momentum 1 = 0. More elaborated calculations for 
massive quarks, more realistic trial wave function and with inclusion of string corrections and 
short-range Coulomb force are in progress now and will be reported elsewhere. 

The parity and charge conjugation of the hybrids are given by 

P = C = (-l) l+s+1 (23) 

for the states with electric gluon (21), and 

P = C = (-l) /+s+1 (24) 

for the states with magnetic gluon (22), where s is the total spin of quarks. So the possible 
quantum numbers for the ground state hybrids are 

jPC = t +! 1 t+ ) 2 ^+, I*-. (25) 

With the chosen trial function these states are degenerate, and we expect this degeneracy to be 
removed mainly by spin-dependent effects. For the string tension o = 0.2GeV 2 the eigenenergy 
of the ground state is estimated to be E = 3.3GeV. 

Now we enter the rather delicate point of setting the absolute scale for the mass of a hybrid. 
It is well-known that in the potential models the large and negative additive constant in the 
potential is needed to fit the meson spectra. The same holds for the lowest qq's in the suggested 
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picture: the negative constant C ~ —0.7 -. — 0.8 GeV is needed to describe the S—,P— and 
D — qq levels. At present there is no reasonable model for this constant term, and the best 
way would be to adjust it by some hybrid candidate mass. In principle, the constant may 
be attributed to the perimeter term in the formula for Wilson loop, or to the hadronic shift 
(creation of additional hadronic loops due to the string breaking). For the qqg system with two 
strings the most naive recipe to account for such mechanisms is that the constant is twice as 
large as for the qq system. Taking this as an educational guess we end up with the mass 



for the ground qqg state. 

The value (26) should not of course be taken seriously. However, it appears to be rather 
close to the mass range where the hybrid candidates are believed to be placed, and to the value 
obtained in the flux tube model [3,4]. Ideologically, our approach and the flux tube model have 
common features: both models account for the transverse motion of the qq string, though the 
idea is realized in different ways. The most important difference is in quantum numbers: there 
is no constituent gluon in the flux tube, so instead of (25) one has [3] 



It is not easy to distinguish between models using the existing data in light quark sector. 
Indeed, at present several hybrid candidates in the mass range 1.5-2 GeV are under discussion: 
1~ + exotic signal from BNL experiment [16], the 0~ + and 2~ + states seen by VES collaboration 
[17], the 7r(1775) state [18] which may be 1 h , but 0~ + or 2~ + assignement is not excluded, 
and additional vector states seen in e + e~ annihilation around 1.5-1.7 GeV [19]. The hybrid 
assignement is attached to these states basing on the very clear signature for hybrid decays 
in the flux tube model [20]: the decay of a hybrid into two S-wave mesons is suppressed, 
and it is really the case for the listed states. Nevertheless, just the same signature exists 
for a hybrid with electric constituent gluon [21,22], that means that the P-odd ground state 
hybrids (25) exhibit the same decay pattern. So to tell one model from another one is to study 
the resonance activity in 0++, l ++ ,2++, l+-/0 +- , l + ",2+-, 1++ channels around 1.5-2.0 GeV, 
a rather challenging experimental task. 

More promising is to use another selection rule, valid only for the decay of constituent 
hybrids [21,22]. There are two decay mechanisms. One proceeds via string breaking, leads to 
the qq — qqg final state and is forbidden for the ground state because of phase space. Another 
is due to the conversion of a gluon into qq pair. The constituent gluon behaves as the massive 
particle, transferring its spin to the created pair, and it means that not only the total angular 
momentum but also the total spin of constituents conserved, giving rise to the set of selection 
rules for hybrid decays, similar to those obtained in [22]. 

To conclude, we have demonstrated that the model for QCD quark-ant iquark string can 
be extended to account for hybrid meson excitations. The arising model of constituent glue 
appears to be compatible with the data on light quark meson spectroscopy. 
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M = 1.7 - 1.8GeV 
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Figure caption 

Wilson loop configuration corresponding to the propagation of the hybrid state. 
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